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ABSTRACT 

We present high spatial resolution maps of ro-vibrational molecular hydro- 
gen emission from the environment of the GG Tau A binary component in the 
GG Tau quadruple system. The H 2 v = 1-0 S(l) emission is spatially resolved and 
encompasses the inner binary, with emission detected at locations that should be 
dynamically cleared on several hundred-year timescales. Extensions of H2 gas 
emission are seen to ~100 AU distances from the central stars. The v = 2-1 S(l) 
emission at 2.24 //m is also detected at ~30 AU from the central stars, with 
a line ratio of 0.05±0.01 with respect to the v = 1-0 S(l) emission. Assuming 
gas in LTE, this ratio corresponds to an emission environment at ~1700 K . 
We estimate that this temperature is too high for quiescent gas heated by X-ray 
or UV emission from the central stars. Surprisingly, we find that the brightest 
region of H 2 emission arises from a spatial location that is exactly coincident 
with a recently revealed dust "streamer" which seems to be transferring material 
from the outer circumbinary ring around GG Tau A into the inner region. As a 
result, we identify a new excitation mechanism for ro-vibrational H 2 stimulation 
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in the environment of young stars. The H2 in the GG Tau A system appears 
to be stimulated by mass accretion infall as material in the circumbinary ring 
accretes onto the system to replenish the inner circumstellar disks. We postulate 
that H 2 stimulated by accretion infall could be present in other systems, particu- 
larly binaries and "transition disk" systems which have dust cleared gaps in their 
circumstellar environments. 

Subject headings: physical data and processes: accretion, accretion disks — stars: 
pre-main sequence — stars: winds, outflows, accretion — stars: formation - 
planetary systems: protoplanetary disks — stars: individual (GG Tau A) 



Introduction 



Planet searches and corresponding imaging follow-up suggest that as many as 25% of 



extrasolar planets reside in binary or higher order multiple stellar systems fe.g.lPatience et al. 



2002 



Eggenberger et al.ll2004l ; iMugrauer et al. I l2007t iTakeda et al.ll2008l ; IMugrauer fc Neuhauserl 



20091 ). This multiplicity fraction is likely a lower limit, since most searches have preferentially 
selected the potential hosts to be single star systems. Such a high multiplicity fraction for 
exoplanet host systems indicates that despite the increased dynamical complexity of young 
multi-stellar systems, which lead to the truncation of individual circumstellar disks and the 
rapid clearing (<1000 yr) of the inner regions of circumbinary(system) disks, reservoirs of 
material frequently survive long enough for planets to form. Stars commonly form in binary 
systems, so the fact that planets exist in th ese systems adds tremendously to the diver- 



sity and prevalence of planets in the Galaxy (Ipuauennov fc Mayor 



Reipurth fc Zinneckerlll993l ; iLeinert et al.lll993 



Simon fc Prato 



1991 



19951 : 



Ghez et a 



Duchene et al 



1993 



20071 ) 



Gas flows through dust-cleared gaps in the disks of binary star systems can also provide 
important clues about the processes that could occur in the embedded planetary systems 
that open gaps. Planetary systems can be considered to be binaries of extreme mass ra- 
tios. Most of the mass gain of gas giant plan ets may occur by means of mass flow through 



the gap surrounding the orbit of the planet fjLubow et al.l Il999l ; iLubow fc D'Angelol 12006 



Ayliffe fc Batdl2010f ) The gas flow through gaps in young planetary systems and onto circum- 
planetary disks plays a crucial role in the formation of the young planets. The physics of 
this process can be better understood with direct observational constraints placed on the 
similar processes that occur in young binary star systems. 



In order to understand the formation of planets in young binary and multiple systems, 
we first must understand the ways in which the stellar cores interact with the gas and 
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dust in their circumstellar environment. Clues to the dynamics of these interactions and 
the structures in the circumstellar /system environment are revealed by spectral energy dis- 
tribution measurements, high-spatial resolution observations, and hydrodynamic simulations 



d Jensen fc Mathieul 19971 



19961: 



Giinther fc Kiev 



20021; 



Pascucci et al. 



Artvmowicz fc Lubowlll994l; iDutrey et al.lll994l ; lArtymowicz Sz Lubow 



20081 ). For instance, simulations showing how 



young binary systems quickly clear gaps in the inner regions of circumbinary disks also de 
scribe how material will fall through these ga ps in "high- velocity streams," eventually crash 
ing i nto the individual circumstellar disks ( lArtymowicz fc Lubowl Il996l ; iGiinther fc Kiev 
20021 ). High-spatial resolution millimeter and optical/near-infrared observations of nearby 
young binary systems like GG Tau A clearly sho w that optically-thin regions exist between 



the circumbinary disk and the stellar components (IDutrey et al 



Roddier et al. 1996: McCabe et al. 2002; Krist et al. 



2002 



2005 



1994 



Guilloteau et al. 1999 



Pietu et al.ll201ll ). However 



these observations also strongly hint at the prese nce of streamers o f material falling through 
the optically-thin, dynamically unstable region ( jPietu et al .1120 111 ). More recently, a near- 
infrared high-resolution image of scattered light from SR 24, a binary system in nearby the 
Ophiuchus star forming region, clearly indicates the presence of circumstellar disks orbiting 
both the primary and secondary and a spiral arm-like str eam of material exte nding some 
1600 AU from central sources toward a circumbinary disk (IMayama et al.ll2010l ). 



The replenishment of the tidally-truncated, inner circumsie/Zar disks with material from 
the extended circumbinary (system) disks is likely to aid the formation of planets in many 
binary and multiple systems. Moreover, the recent discovery of Kepler-16b, an exoplanet that 
lies within the same orbital plane as its c entral binary host , demonstrates the potential for 
planets to form within circumbinary disks (IDoyle et al.ll201ll ). High-spatial resolution images 
capable of resolving emission from structures within a circumbinary region and near the inner 
circumstellar disks will be essential for studying the planet-building process in such systems. 
Constraining existing hydrodynamical models, which form our basis for understanding the 
evolution of young multiple systems and their host planets, would be a meaningful step in 
unraveling the complexity of star and planet formation in the binary environment. 

As an exciting step toward this goal, we present the first spatially resolved image of 
molecular hydrogen gas located within the optically-thin, unstable region of the GG Tau A 
system. This image provides strong support for the existence of suspected accretion "stream- 
ers" in this system. The material falling onto the inner ambient gas appears to create a shock 
front located in the dynamically unstable circumstellar environment of the stars. We compare 
the ratios of multiple H 2 lines to those predicted by models of shocked molecular hydrogen 
gas and conclude that the emission is consistent with shock heating. We use these data to 
estimate the mass and physical conditions of the emitting gas residing within the cleared 
region of the GG-Tau A's circumbinary disk. Finally, we discuss the phenomenon of shocked 
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infalling gas in a binary system and the implications for similar systems such as transitional 
disks that potentially harbor planetary companions in their inner dust-cleared gaps. 



2. Observations 

Observations were obtained using the Near IR Integral Field Spectrograph (NIFS) at 
the Gemini North Frederick C. Gillette Telescope on Mauna Kea, Hawaii. NIFS is an image 
slicing IFU fed by Gemini's Near IR adaptive optics (AO) system, Altair, that is used to 
obtain 3-dimensional x, y, wavelength (A) spectral cubes at spatial resolutions o f <0'.'l with 



a 2 pixel spectral resolving power of R~5300 at 2.2 /im ([McGregor et all 120031 ) . The NIFS 
spatial field is 3" x3" in size and the individual IFU pixels provide O'.'l x0."04 spatial sampling. 
Data were obtained at the standard X-band wavelength setting for a spectral range of 2.010- 
2.455 jum. GG Tau A was bright enough to serve as its own wavefront reference stars for 
the facility adaptive optics system to guide at a 1000 Hz rate. All observations used for the 
final co-added datacube were acquired in natural seeing of better than ~0?85 for excellent 
AO correction. 

The data for this project was acquired in queue observations during Fall 2009 under pro- 
gram ID GN-2009B-Q-40 through thin cirrus and good seeing. Observations of GG Tau A 
were executed over four different nights. The data reduction and analysis of th e IFU dat- 



acub e s on G G Tau was identical to our past integral field spectroscopic studies ( iBeck et al. 



20081 . 120101 ) . The individual reduced GG Tau A datacubes were re-sampled to 0'.'04 x 0'.'04 
size pixels in the spatial dimensions. The exposures each had 40 second integration times, 
and the data was acquired with small dithers that were the same size as the width of the 
NIFS IFU slicing mirrors. The purpose of dithering the stars at integer slice mirror widths 
was to sample the stellar point spread function (PSF) in the cross dispersion direction in dif- 
ferent pixels to improve the spatial sampling and mitigate the effect of the NIFS rectangular 
spatial elements. The final datacube on GG Tau A was constructed by shifting-and-adding 
78 individual dithered exposures for a total on-source exposure time of 3210 seconds. To 
verify the i^-band flux of the GG Tau A system and calibrate the imaging spectroscopy, sup- 
porting photomet ry was acquired in S ep. 2009 using the WIYN High-Resolution Infrared 



Camera (WHIRC; iMeixner et al.ll2010l ). using the K s filter and the United Kingdom Infrared 
Telescope (UKIRT) faint standard FS 18 for flux calibration. The derived K s magnitude 
for the combined GG Tau A brightness is 7.7±0.1, which is consistent with all published 
photometry of the system. 



By inspecting the wavelength calibration and sky emission line positions, we estimate 
that the absolute velocity accuracy of our individual IFU cubes for a given night of obser- 
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vations is ~9-12 km s _1 (roughly one third of a spectral pixel, or one sixth of a resolution 
element). For our observations, the accuracy of the NIFS IFU grating wheel repositioning 
was repeatable only to ~ ±0.75 pixels over the four nights of our observations. Hence, 
because of very weak signal from H 2 emission in GG Tau A, data acquired over multiple 
nights was not subsampled and shifted and co-added in the wavelength dimension to pre- 
serve velocity resolution. As such, we estimate that our final emission line velocity centroid 
accuracy is affected by the grating repositioning and the straight co-addition of data that is 
slightly shifted in wavelength. As a result, the estimated absolute velocity accuracy of the 
GG Tau A data is at the level of 0.7 to 1.0 pixels, or 20-28 km s" 1 . 



3. Results 

Figure 1 presents the image of 2.12 /xm v = 1-0 S(l) H2 emission in the environment of 
GG Tau A. This image was constructed by summing six spectral pixels centered on the H 2 
emission at each spatial location in the field, after the continuum level was subtracted using 
a linear fit to the stellar continuum flux on each side of the emission feature. GG Tau Aa, 
the southern star, is centered in the image. Overplotted in blue are the 2/im continuum 
contours for both sources. The H 2 emission is distributed around the binary, with obvious 
extensions of the emission to ~70 AU (0'.'5) distances from GG Tau Aa. The strongest H 2 
emission arises from north-east of the system, about 0'.'2 away from the peak continuum 
position of GG Tau Ab. Given the projected separation of the GG Tau A binary, the H 2 
gas emission arises from dynamically unstable locations near the stars. In the image, the 
apparent decrease in H 2 emission strength at the positions of the stars is likely not genuine, 
it results from decreased sensitivity to the faint H 2 over the very bright stellar flux. If H 2 
at the level of the peak emission seen in Figure 1 existed at the location of GG Tau Aa, it 
would be less than 0.5% of the continuum flux and would be extremely difficult to detect 
after the continuum subtraction process. 

A binned 2x2 image of the H 2 emission from the environment of GG Tau A is shown in 
Figure 2a. In this image, the brightest H 2 is saturated to emphasize the low level emission 
in regions surrounding the young binary. H 2 "streamers" of emission are seen extending 
out to ~0'/9 distances from the central stars. Three positions are highlighted in the image - 
Positions A, B and C. To the right in Figure 2b, the 1-D spatially integrated spectra (from 
0'.'12 radius apertures) of the 2.12 /xm v = 1-0 S(l) line are presented. The integrated signal- 
to-noise ratios are ~24, 4 and 16 for the A, B and C spectra, respectively. The velocity of the 
H 2 emission in the A and B spectra are centered at ~-10 km s _1 . Given the ~ ±20-28 km s _1 
absolute accuracy of our velocity calibration, this is not a statistically significant difference 
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from the systemic radial velocity of GG Tau A of 18 km s 1 f lHartmann et al.lll986f ). These 



results are consistent with the previous high spectral re solution (R~60,000) detection of H 2 



emission from the GG Tau A system flBary et all 120031 ) . The spectrum extracted at the C 



position has a double-peaked profile shape, with a significant redshifted contribution to the 
line flux. 

Figure 3 presents a spectrum of the H 2 emission summed from spatial locations that 
have greater than 4.0xl0~ 17 W m~ 2 line flux in the v = 1-0 S(l) line. In addition to the 
v = 1-0 S(l) line at 2.12 /mi, five other H 2 emission transitions are detected, including the 
v = 2-1 S(l) line at 2.24 /mi. The measured value of the v = 2-1 S(l) / v = 1-0 S(l) line ratio 
is commonly used to determine the physical conditions of the emitting H 2 gas. Assuming 
that the gas is in LTE, the Boltzmann relation and line ratio provide an estimate of the 
temperature of the emitting gas. For the spectra in spatial locations nearby to GG Tau A 
(Figure 3), we measure a value of ~0.05±0.01 for this line ratio and determine a correspond- 
ing gas temperature of ~ 1700 ±100 K. At this temperature and average ~40 AU distance 
from the central stars, it is unlikely that this gas is heated by blackbody emission from the 
stars. Given that the line ratio is also consistent with shock heating or high-energy photo- 
excitation in a dense medium (n# > 10 4 cm -3 ), both processes are potential mechanisms 
for stimulating the observed emission. However, given the spatial extent of the emission, the 
lack of detectable X-ray emission and relatively weak UV flux from the binary, we conclude 
that shock heating is the most likely excitation mechanism for the H 2 . 

The detection of v — 2-1 S(l) emission at ~0'.'2 distances from the central stars is a very 
high-contrast measurement. The peak flux in the brightest region of v = 1-0 S(l) emission is 
0.6% of the continuum flux at the position of GG Tau Ab. The detected v = 2-1 S(l) emission 
is 5% of the flux of the 1-0 S(l) line. This discrete line emission is detected at 5a confidence, 
with a contrast ratio of 3xl0 -4 with respect to the continuum flux level of GG Tau Ab. Our 
total emission line sensitivity is better than this because the stellar PSF of only the closer 
GG Tau Ab star in the binary is considered. We estimate the overall sensitivity at the ~0'.'2 
distance from both stellar PSFs to be ?«8xl0 _5 . At this sensitivity level, we are successfully 
probing the parameter space of the next generation planet survey instruments, such as 
those using extreme c oronographic adaptive optics and non-redundant masking techniques 



( jBeichman et al.ll2010f ). This is very difficult, and is made possible primarily because we are 
searching for discrete, narrow emission lines where the adaptive optics halo structure can 
be accurately subtracted from the stellar point spread function over the small wavelength 
range. Our detection of v = 2-1 S(l) emission in the GG Tau A system highlights the success 
of using adaptive optics fed IFU spectroscopy to spatially resolve emission line morphologies 
within the inner 30 AU circumstellar environments of bright young stars. 
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Discussion 



The GG Tau system was first spatially resolved into a hierarchical quadruple by lLeinert et al 

( 1199 lh . At the time of those observations, the stars in the GG Tau A component of the 
quadruple had a spatial separation of 0'.'255 at a position angle of +9°. Our observations, 
acquired ~20 years later, show a system separation of 0.251±0.008 at position angle 334±1 
(measured east of north). In the twenty years since its discovery as a binary, the GG Tau A 
stars have cha nged their or ientation by nearly 35 degrees, just under one tenth of the total 



angular orbit ( )Kohlerll201ll ). The fact that the separation has remained largely constant over 
a wide arc of the orbit means that the system is viewed close to pole on and either the orbit is 
nearly circ ular, or the c omponents are near periastron where the apparent separation would 



at 



White et al. 


1999; 


Dutrev et al. 


1994 



Guilloteau et al.lll999l ). lHartmann et all ( 119861 ) found the combined GG Tau A system to 



have a vsini of less than 10 km s _1 , implying that the stellar components may also be viewed 
close to pole-on. Though, in order to explain the large ~180 AU extent of the inner edge of 
the circumbinary ring, dynamical modeling of the system implies that the orbital eccentricity 
of GG Tau A may lie in the range of 0.2-0.4 and the stellar bina ry orbit may not be coplanar 



with the circumbinary ring (IBeust &: Dutreyll2005t lKohlerll201lh . Yet, when the full GG Tau 



system is modeled as a quadr uple then the need for non-coplanarity of GG Tau A with the 



circumbinary ring is removed (IBeust &: Dutreyl 120061 ) . 



The peak of the detected molecular hydrogen emission in the GG Tau A system is located 
at a position that is ~0'. / 2-0 / .'25 away from both GG Tau Aa and Ab. Given the separation 
of the components of the binary of 0'.'25, the tidally truncated inner circumstellar disks will 
not extend to this distance from the stellar positions. Hence, H 2 detected at this position 
is too extended to be associated with either of the circumstellar disks that exist around 
the stars. Hydrodynamic models simulating the clearing of the inner regions around young 
star binaries clearly suggest that material at the observed lo cation of H? emission should 



be dynamically cleare d on timescales of a few hundred years (lArtymowicz fc Lubowlll996 



Giinther fc Kleyi 120021 ) . Though, the presence of material within the dynamically unstable 
region is not unexpected, given that GG Tau system possesses several indicators of accretion 
activity and the inner disks must be replenished to persist. 

The circumbina ry ring encircling GG Tau A was first detected in scattered light by 
Roddier et al.l ( 119961 ). Since the time of the first scattered light images, higher sensitivity 
observations have revealed an inner distribution of dust within the circumbi nary ring of 
GG Tau A, extending to ~0'.'5 north of the inner binary ( iDuchene et al.l 120041 ) . This inner 
dust distribution corresponds in spatial location to the northern extension of H 2 emission 
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that we see. Our detected extended H2 emission can not merely be scattered line flux from a 
compact, centrally excited H 2 emission environment because a central emission region would 
need to be hundreds to thousands of times stronger than our detected extended line emission 
region. To produce our measured extended H 2 line flux, a central emission region would need 
to be at a level of more than a 30% of the continuum, and would be easily detectable over 
the central stellar flux. Moreover, if our detected H2 line emission arose from scattered line 
emission, we would also expect to detect H 2 at the circumbinary ring, which is the location 
of the strongest and most easily detected scattered flux. We do not detect appreciable H 2 
emission from the circumbinary ring, to a detection limit of 8.0xl0~ 19 W/m 2 /pixel. The 
brightest region of H 2 emission we detect seems to be spatially coincident with the inner 
region of dust-scattered light revealed in broad-band imaging, though the H 2 emission must 
be excited locally within this environment rather than scattered. 

Various high contrast observations of the GG Tau A system have reported the detection 
of "spokes" or "streamers" of material extending from the inner edge of the circumbinary 
ring inward toward the central binary system. These observations have been met with some 
level of skepticism because while ground-based studies seem to show some weak evidence for 



the streamers, imag es from the Hubble Sp ace Telescop e typic ally did not ( jKrist et al.l 12005 
McCabe et al.ll2002l ). For the first time, iPietu et al.l ( 1201 ll ) detected a possible streamer 



of dust emission in the mm wavelength map of the circumbinary ring around GG Tau A. 
This streamer of dust was detected at a ~2cr level of confidence, and it seems to extend 
from the northern edge of the circumbinary ring inward to the north side of the central 
binary. Figure 4 presents an image of our H 2 ga s emission a r ound the GG Tau A system 
with the 1.1 mm dust continuum contours from IPietu et al.l (120111 ) overplotted. Our H 2 
map also shows evidence of "streamers" of emission, and several of the H 2 extensions are 
detected at high signal-to-noise. Interestingly, the location of the inner edge of the mm 
dust streamer corresponds precisely to the location where we detect the peak emission in 
molecular hydrogen gas. The presence of the peak H 2 emission at a coincident location 
with the inner dust strea mer lends further support th at the streamer in the dust map is a 
real feature. Additionally, iGuilloteau fc Dutreyl (120001 ) report the existence of CO emission 
from these regions, which would also support the existence of streaming gas interior to the 
circumbinary dust ring. 

In the environments of young stars, most excitation mechanisms considered for ro- 
vibrational H 2 excitation involve stimulation by X-ray o r UV fl ux from the central star, or 
shock exc itation by mat e rial in an outflow. iBary et al.l (120031 ) using the X-ray excitation 
models of iMaloney et al.l (119961 ). ignored the complexities of the binary system and assumed 
a centrally located X- ray/UV source incident on material located in a standard disk model 
(IGlassgold et al.ll2000l ) to predict a range of plausible H 2 line fluxes for the GG Tau A system. 
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Using log L x = 29 A ergs s~^_, Bary et al. (2003) assumed that the X-ray energy was deposited 
in the upper atmosphere of a singular disk between 10 and 30 AU from the X-ray source. In 
the best case scenario, the X-rays were capable of sti mulating rou g hly 70 % of the observed H 2 
line flux in a high-resolution spectrum of GG Tau A. iBary et al.l ( 120031 ) further assumed the 
UV flux from GG Tau to be similar to that of TW Hya and estimated the contribution of UV 
fluorescent H2 emission to be substantially less than the X-rays. However, they concluded 
that the combination of the limit of X-ray ionization and UV fluorescence might account for 
a significant portion of the H 2 line flux in the GG Tau A system. With the spatial resolution 
and wavelength coverage of NIFS, we find that the gas is not distrib uted in the manne r 
predicted by high -energy photon disk excitation models described by IBary et al.l ( 120031 ); 
Bary et al.l (120081 ). The H2 emitting gas extends farther from the sources than predicted 
and appears to display a morphology better described by gas participating in an accretion 
streamer than the central compact circumstellar disks illuminated by X-rays. In addition, the 
v = 2-1 S(l) / v =1-0 S(l) line ratio of 0.05±0.01, measured at ~40 AU distances corresponds 
to T ~1700±100 K. This ratio can be accurately modeled by C-type shocks with reasonable 
pre-shock gas densities and shock velocities (see below; iLe Bourlot et al.l 120021 ) . In light 
of these considerations, the lack of detectable X-ray flux from the GG Tau A system, and 
the spatial coincidence of the northern H2 gas with the possible dust stream (Figure 4), 
we conclude that this gas is most likely being shock excited by infalling material crashing 
into gas distributed throughout the circumstellar domain than stimulated by high-energy 
photons. 

The measured emission line flux of the v = 1-0 S(l) transition an d the estimated gas 
excitation temperature coupled with the partition function analysis from lHerbst et al.l (119961 ) 
allow us to determine the molecular hydrogen column density: 



N tot = 4.55 x 10 16 • T ■ (e 600 °/ T - l)e 956 / T • h ■ 10 a4A > 



where I\ is the measured line intensity (in erg cm 2 s 1 sr x ), T is the gas temperature 
(assuming LTE), a nd A^ is the extinction to the emitting region, which we adopt to be 0.7 



( [White et al.lll999l ). The derived column density of H 2 in the brightest knot of emission is 
N tot = 1.6 x 10 15 cm -2 . The H 2 emitting region sampled here was a square pixel with a 
O'.'l size. Hence these dimensions and the calculated N tot implies a mass of ~ 10 -12 M in 
molecular hydrogen in this aperture, assuming a 140 pc distance to the Taurus molecular 
cloud. Assuming comparable column densities in regions of the H 2 circumbinary material 



1 Bary et al.l ([20031 ) did not note that this value was taken as an upper limit since GG Tau A was not 
detected by ROSAT 
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encompassing the system out to ~0'.'3 distances, we estimate that the mass in H2 around 
GG Tau A is at least several xlO -11 M . Using the stellar paramete rs from GG Tau Aa an d 



Ab from IWhite et al.l (119991 ) with the accretion rate relations from iMuzerolle et al.l ( 119981 ). 
we have used our measured Br7 line luminosity to derive mass accretion rates onto the stars. 
These accretion rates are 2.7xl0~ 8 M yr" 1 and 2.4xl0" 8 M Q yr" 1 for GG Tau Aa and Ab 



respe ctively, consistent with past mass accretion rate estimates for the system (IWhite et al. 



19991 ). Not only does dynamical theory predict that the H 2 should not be stable at our 
detected locations, but the mass accretion rate onto the stars is several orders of magnitude 
larger than the estimated mass of emitting H2 in the circumbinary environment. At this level 
of mass accretion, the extended H 2 material must be continually replenished or it would be 
accreted rapidly onto the central disks and stars. 

We have compared the detected fluxes of the v = 1-0 S(0), S(l), S(2) an d v = 2-1 S(T) H 2 



l ines f rom the spectra of GG Tau (Figure 3) with the C-type shock models of lLe Bourlot et al. 



(120021 ). We probed a range of pre-shock medium densities from 10 3 to 10 7 cm -3 and shock 
velocities from 10 km s" 1 to 70 km s -1 . The very slow shocks (<10 km _1 ) are ruled out 
with confidence, at these low velocities the shocks cannot produce the measured excitation 
temperatures greater than 1000 K. Extremely fast high-density shocks also cannot reproduce 
the measured line fluxes, as these models result in greater excitation temperatures measured 
by the v = 2-1 S(l) / v = 1-0 S(l) line ratio. We find that the optimal solutions in the 
models are for a pre-shock medium in the density range of 10 3 — 10 5 cm -3 with a shock 
velocity of 20-50 km s _1 . In GG Tau A, the accretion infall shock velocity would result from 
a combination of the gas free-fall velocity onto the system (~10 km s _1 ) with a component 
from the Keplerian motion of material in the inner system. With these constraints, only the 
model calculations of shock velocities in the lower 20-30 km s _1 range would be consistent 
with the expected gas accretion kinematics for GG Tau A. 

We speculate that H 2 can also be stimulated into emission by mass accretion infall in 
other systems, such as young star binaries and systems with "transitional disks" that have 
dust cleared disk gaps that are thought to trace the onset of planet formation. In these latter 
systems, gas and dust from the outer distribution of disk material could accrete onto the 
inner regions, and H 2 might be stimulated into emission from the shock excited free-falling 
material, which could have velocities in the 20-40 km s _1 range. For example, the GM Aur 
transition disk system has a dust-cleare d gap in the inner d i sk that extends out to a distance 
of ~20-25 AU from the central star (jDutrey et al.l 120081 ; iHughes et al.l 120091 ). Accretion 
streams of gas that fr ee-fall from the ~20 AU outer ring onto an optically thin inner disk 
in the <1 AU region ( jEspaillat et al.ll2010l ) could have a velocity at the 30-40 km s _1 level, 
which would be fast enough to shock-excite the H 2 into emission. 
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In addition to the peak of H2 emission coincident with the inner dust streamer, ex- 
tensions of H 2 emission exist out to ~0'.'7 (100 AU) separations from the central stars in 
GG Tau A (Position B in Figure 2). The detected H 2 gas itself seems to exhibit "streamers" 
of emission that radiate outward from the central system. More sensitive data at both mm 
and near-infrared wavelengths could clarify the correlation between the streamers and the 
dust and gas species present within the system. Further study of the spatially resolved ma- 
terial apparent in the GG Tau A ring system would lead to a more detailed understanding 
of the general physics of mass accretion onto a central stellar binary. While this knowledge 
is directly applicable to other young star binaries, it also can serve to constrain the Kep- 
lerian motion in a broader base of models such as those that describe mass accretion onto 
proto-planets within circumstellar disks, as well as binary black holes in galaxy mergers and 
binary stellar remnants. 
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Fig. 1. — A Continuum subtracted grey scale image of H2 emission in the environment of GG 
Tau A binary, with continuum contours at 0.1, 0.4 and 0.7 times the peak continuum flux 
of GG Tau Aa overplotted in blue. The 0.7 and 0.4 contour levels were selected to highlight 
the continuum flux positions of the GG Tau Aa and Ab stars in this young binary. The H2 
emission is distributed around the binary, with extensions to the emission out to ~100 AU 
(0."7) distances. The apparent decrease in H 2 emission at the positions of the stars is likely 
not genuine, it results from uncertainty in the continuum subtraction process to detect the 
faint H2 over the bright stellar flux. 
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Fig. 2. — To emphasize the detection of low level line emission, a binned 2x2 image of the H2 
emission from the environment of GG Tau A is presented in a). In this image the display is 
saturated to show the emission surrounding the young star binary. Three positions - Position 
A, B and C - are highlighted in color in the image. To the right in b), the 1-D spatially 
integrated spectra of the 2.12/iin v = 1-0 S(l) line at the A, B and C positions are presented. 
The H 2 line emission is detected at significant signal-to-noise at spatial positions of more 
than 100AU from the central stars. 
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Fig. 3. — A spectrum of the H2 emission summed from regions of flux that are brighter 
than 4.0xl0 -17 W m -2 (as seen in Figure 1). In addition to the v = 1-0 S(l) line at 
2.12/im, five other H 2 emission transitions are detected, including the v = 2-1 S(l) line at 
2.24/im. The v = 2-1 S(l)/v = 1-0 S(l) ratio we detect from the gas around GG Tau A is 
~0.05±0.01, implying the presence of warm, T ex ~1700K gas at extended ~40 AU distances 
from the system. Residual stellar photospheric absorption features of Al+Mg (2.11/im), Na 
(2.206/2.208^m), Ca (2.26 fan) and CO (2.29 - 2.38/im) are seen in this spectrum. 
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Arc Seconds 

Fig. 4. — The continuum subtracted H2 image in the environment of GG Tau A (grey) with 
overplotted contours of the 1.1mm dust continuum emission (blue) and the 2/iin continuum 
flux (red). The brightest peak of molecular hydrogen emission corresponds to the spatial 
location of the dust "streamer" that is seen to extend within the circumbinary ring, which 
implies that both features are real and the H2 may be shock excited by mass accretion infall 
onto the binary stars. 



